Force fluctuation in a driven elastic chain 
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We study the dynamics of an elastic chain driven on a disordered substrate and analyze numer- 
ically the statistics of force fluctuations at the depinning transition. The probability distribution 
function of the amplitude of the slip events for small velocities is a power law with an exponent 
r depending on the driving velocity. This result is in qualitative agreement with experimental 
measurements performed on sliding elastic surfaces with macroscopic asperities. We explore the 
properties of the depinning transition as a function of the driving mode (i.e. constant force or con- 
stant velocity) and compute the force- velocity diagram using finite size scaling methods. The scaling 
exponents are in excellent agreement with the values expected in interface models and, contrary to 
previous studies, we found no difference in the exponents for periodic and disordered chains. 
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I. INTRODUCTION 

The dissipative motion of an elastic line in a random 
potential is an interesting example of a nonequilbrium 
interacting system and is relevant for several phenom- 
ena in condensed matter physics. Examples are numer- 
ous and include the motion of magnetic interfaces in fer- 
romagnetic materials solid friction |||-|^, wetting 
charge density waves fluids in porous media 
1^1, vortex dynamics in high temperature superconduc- 
tors ]lO|-p^, cracks and dislocations. These systems are 
characterized by a dynamic phase transition ruled by the 
interplay between quenched disorder and elastic interac- 
tions. 

Due to the effect of the disorder, an elastic chain at 
zero temperature is pinned when the applied force is be- 
low a critical value Fc'. after a sufficiently long time, in- 
dependently of the initial conditions, the chain reaches a 
configuration where no movement is possible. For F > Fc 
the chain can escape from any pinning configuration and 
moves with constant average velocity. When F is close 
to Fc the motion is dominated by collective effects and 
the depinning of a single bead produces a large reorga- 
nization of the chain. In other words, for F ^ Fc the 
system is critical and the motion of the beads is highly 
correlated. 

An elastic chain moving in a disordered potential is a 
useful model to understand some general features of slid- 
ing friction in particular of the experiment reported in 
Rcf. []l3|-^ , done using two artificial surfaces with con- 
trolled roughness and elasticity. Beads of diameter 2mm 
were randomly put inside an elastic matrix, with a max- 
imum roughness of 0.5mm. The two surfaces were then 
displaced against each other at constant velocity and the 
friction force was measured, varying the elasticity of the 
matrix and the driving velocity. The distribution of the 
amplitude of the slip events is generally found to decay as 
a power law at small velocities, suggesting the presence 
of an underlying critical point. The exponents charac- 
terizing the power law distribution are found to decrease 



with the applied velocity, in analogy with other driven 
systems such as domain walls in ferromagnet . 

Several variant of the chain model can be studied in or- 
der to reproduce the experiments: periodic or disordered 
H arrangements of beads, on a rigid or elastic 
substrate. Simulations of a one dimensional, periodic or 
disordered, chain over a rigid disordered potential have 
been performed by Cule and Hwa jj]. The measurement 
of the velocity and the roughness exponents seem to in- 
dicate that periodic and disordered chains are described 
by two different universality classes. These and other 
simulations are performed considering a constant applied 
force, while experiments are performed driving the sys- 
tem at constant velocity. 

Here we consider explicitly the second case and mea- 
sure the force fluctuations as a function of the applied 
velocity. The distribution of the slips events is power law 
distributed, and characterized by an exponent t, which 
appears to decrease with the applied velocity, in agree- 
ment with the experiments. In the limit of low velocity 
the exponent r can be related to the critical exponent 
obtained tuning the applied force as first discussed in 
Ref (see also ) . We measure the fluctua- 

tion in the position of the beads in the constant force and 
constant velocity cases and obtain the same roughness ex- 
ponent C. Next we evaluate the force velocity diagram, 
using finite size scaling to locate the critical force and 
compute the exponent /?. The values of the exponents 
are consistent with scaling relation and in good numerical 
agreement with the exponents of interface depinning, but 
disagree with previous simulations for a periodic chain in 
a random potential [Q . In order to confirm this conclu- 
sion, we simulate the motion of a disordered chain and 
see no evidence for the existence of two different uni- 
versality classes for periodic and disordered chains, in 
disagreement with the conclusions of Ref. Q . 

The paper is organized as follows in Sec. |l^ we intro- 
duce the model, and in Sec III we define the critical ex- 
ponents and discuss some scaling relations. In Sec. IV A 
we present the numerical results obtained at constant 
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velocity, in Sec. |IVB we discuss the constant force case 
studying the scahng behavior close to the depinning tran- 
sition. In Sec. we summarize the main results of the 
paper. 



II. MODEL 



We consider the overdamped dynamics of a one dimen- 
sional line of elastically coupled beads, driven on a disor- 
dered substrate. The disordered potential is a succession 
of identical Gaussian potentials, randomly distributed in 
space. The beads can be driven directly, applying a con- 
stant force, or indirectly coupling them to an interme- 
diate spring which is pulled at constant velocity. The 
equation of motion is. 



+ n^i) + f{uit)) + F, (1) 



where ri{t) is the position of the bead i at time t, rj is the 
coefficient of viscosity, D is the stiffness of the elastic line, 
F is the driving force, and f{x) is a random force, due to 
the contribution of the ensemble of pinning centers. We 
model the random force by the sum of N derivatives of 
Gaussian potentials located on the pinning site 
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where C represents the strength of the disorder, a quan- 
tifies the width of the wells, and is the location of the 
pinning site i, which we chose to be Poisson distributed. 
The equation of motion is integrated numerically using a 
fourth-order Runge-Kutta method. 

The interplay between disorder and elastic interactions 
in our model can be understood computing the Larkin 
length II ||2^]. For distances smaller than Z^, the beads 
are interacting strongly and the chain moves coherently, 
while for distances larger than 1^ the random forces be- 
come dominant and the chain deforms considerably. The 
Larkin length can be estimated considering both the ef- 
fect of the rigidity of the line and the strength of the 
disorder and for our model it is given by m 
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where a is the distance between the beads and p is the 
density of pinning centers. In order to analyze the critical 
properties of the system, we have to consider the limit 
where the Larkin length is larger than the mean distance 
between the beads and the dynamics is governed by the 
collective motion of the beads. To this end we carefully 
choose the parameters of the model so that II a. 



III. SCALING RELATIONS 

Depending on the method used to drive the chain the 
measured quantities change, but the corresponding criti- 
cal exponents can be related by scaling relations. Here we 
summarize the scaling properties of the depinning tran- 
sition in the case where the beads are driven by springs 
of stiffness K pulled at constant velocity V and in the 
case where they are submitted to a constant force F. 



A. Constant velocity driving 

Friction experiments are generally performed under a 
constant velocity driving. In some cases, a traction ma- 
chine turning at constant velocity is coupled by a spring 
to the sliding system. In other cases, an effective spring 
coupling is due to the elastic deformation of the material 
driven imposing a constant strain rate far from the slid- 
ing interface. To simulate constant velocity driving 
we attach each bead to a spring of stiffness X , so that 
the force is given by 



F = K{n{0) + Vt-r,{t)) 



(4) 



where V is the applied velocity. When the velocity V 
and the stiffness K are small, the motion displays large 
fiuctuations. In particular, in the limit V and 
-ftT — > the system reaches the depinning transition and 
the force fluctuates around Fc- This feature is common 
to other non equilibrium critical phenomena, such as ab- 
sorbing state phase transitions and self-organized criti- 
cality||,g. 

The distribution of the friction force fluctuations is di- 
rectly related to the size distribution of slip events Ax of 
the chain, since A_F = KAx. Here x = and the slip 

is defined by a drop in the measured friction force (See 
Fig. 1). Close to the depinning transition, we expect that 
the distribution of Ax decays as 



P{Ax) - g-'' h{Ax / Axq) 



(5) 



where h{x) is a scaling function and Axq the cutoff 
value. The value of the cutoff depends on various pa- 
rameters, such as the system size employed in the simu- 
lation. Clearly the avalanche size cannot be greater than 
the total length of the line. Furthermore, we expect that 
the stiffness of the springs K and the driving velocity V 
will in general change the value of the cutoff. In the low 
velocity limit and for large enough system sizes, K be- 
comes the dominant parameter that determines the value 
of the cutoff g. 

The scaling of the cutoff with K can be used to evalu- 
ate the roughness exponent C, using the relation 
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This relation can be obtained noting that the chain can 
not jump over a distance larger than , which represents 
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the length for which the elastic interaction term is over- 
came by the restoring force due to the springs: £'Co~^ ~ 
K. Then, using the scaling relation for Axq ^ <^q, we 
obtain Eq (|) [||J|]. This result implies that P{Ax,K) 
should satisfy the scaling form 



P{Ax, K)K-^^'^ = H{t), t EE AxK'^/^ 



(7) 



Eq. (H) is used to compute t and Eq. (^) and Eq. (Q) are 
used to compute the exponent C,. 

The exponent C, can also be evaluated directly using 
the scaling of the fluctuations of the displacements with 
K . Defining the relative displacements of the beads as 
Ui{t) = Tiit) — [Vt + TiiQ)), the fluctuations can be quan- 
tified by 
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where mi(t) = X)i=o ^*(^)/-^' ^^'^ ^ number of 

beads. The roughness W scales with the correlation 
length ^0 as ~ ^Q, and since ^ K^^l"^ 



(9) 



Eq. (|^) and Eq. (|p) have the same origin and can inde- 
pendently be used to estimate C,. 



B. Constant force driving 

When the system is driven at constant force we expect 
a depinning transition as a function of F. For F > Fc, 
the chain moves with constant average velocity v defining 
an exponent P 



{F~F, 



(10) 



Close to the depinning transition the motion is very ir- 
regular, and large regions of the chain move collectively. 
The correlation length diverges at the transition as 



i^{F^Fc 



(11) 



In order to estimate the critical exponents (3 and v we 
employ a particular finite size scaling method |1 1| , in 
analogy with absorbing state phase transitions ]22|]. 

We first compute the critical force analyzing the de- 
cay of the average velocity with time for different system 
sizes. For finite systems, the average velocity reaches a 
quasi steady state v{F, L). When F > Fc we expect that 
as the size L — > c», w(F, L) approaches a non vanish- 
ing value given by Eq. (|lO|), while it decays to zero for 
F < Fc- At the depinning transition we expect that 



Once F(. is known with good precision, we can measure 
directly the exponent j3 from Eq. (|lo|). As in the con- 
stant velocity case, we can evaluate the roughness expo- 
nent, measuring the width at Fc as a function of L, which 
should scale as 



W{L) - L^. 



(13) 



,{Fc,L)^L-P'\ 



(12) 



For a periodic chain Ref. Q reported (i ~ 0.4 and C — 1-5, 
while for a disordered (3 ~ 0.25 and C — 1-2. These 
last value are consistent with interface depinning that in 
d = 1 yields (3 ~ 0.25 and C 1-25 [fl. 



IV. NUMERICAL RESULTS 
A. Constant velocity. 

The primary interest of this study is to compute the 
exponent r which characterizes the collective motion of 
the particles at F ~ Fc, and in particular its dependence 
on the driving velocity. We note that the avalanche ex- 
ponent was found to decrease with the driving velocity 
in the Barkhausen effect, due to the motion of domain 
walls in a ferromagnet The same effect was observed 
in the friction experiments reported in Ref [jisj . 

In order to reach the scaling regime, we progressive! 
decrease V and K and compute the friction force. Fig 
shows two typical plot for the friction force as a function 
of the position of the line, in Fig. 0a the driving velocity 
is 1/ = 0.05 and in Fig. |l|b F = 5. As the driving velocity 
increases the friction force becomes smoother, and in the 
limit F ;2> 1 we obtain a viscous behavior (F ~ rjV) with 
small relative fluctuations. On the contrary, for small 
velocities the dynamics is jerky: the force increases with 
time until the beads are sufficiently stressed so that the 
chain depins decreasing the force. In this case, the fric- 
tion force displays a characteristic stick-slip pattern. 

We thus measure the friction force drops AF, or the 
slip sizes AX = AF/K and analyze their distribution. 
The distributions are averaged over 20 realizations of the 
disorder for V = 0.05 and 100 realizations for V — 10; 
in all cases the system was composed of F = 1000 beads, 
and the disorder was produced by = 20000 pining 
sites, Poissonian distributed. The value of the exponent 
T is obtained by a direct fit of the linear part of the distri- 
bution plotted in a log-log graph. The results are shown 
in Fig. y, the main graph presents the log-log plot of the 
probability distribution function of the jumps for various 

V and in the inset we report the value of r as function of 
the velocity. We see a slow decrease of the exponent when 
the velocity increases. This result is in good qualitative 
agreement with the experiments reported in Ref. [ p^ . 
The value of r for ^ agrees well with the expo- 
nent obtained in elastic line depinning under quasistatic 
conditions |25| (see table I). 

Fig. ^ shows that the cutoff AXq clearly depends on 

V when K is hold fixed. We can quantify this variation 
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and the result is reported in Fig. ||. For small veloci- 
ties {V < V*) the cutoff is a constant {i.e it does not 
depend on V) which in principle depends on while 
for high velocities the cutoff decreases with V, roughly 
as a power law. Next, we study the behavior of the cut- 
off when V < V* as K is varied. In Fig. || we shows 
the distribution of slip sizes for various K for very small 
driving velocity. We see that the cut off increases as K 
is decreased. Fig. || shows the collapse of the curves after 
the rescaling with K, in accordance with Eq. (|^). 

We also measure the roughness exponent of the system 
following Eq. (||) and the result is reported in Fig. ^. We 
obtain with both methods ^ — 1.26, which is consistent 
with the numerical value found in interface depinning 
p6t , for a disordered chain and with a recent two 
loop renormalization group calculation p^ ]. This value 
nevertheless disagrees with previous results on a periodic 
chain 

The simulations of Ref. Q for periodic and disordered 
chain, suggest the presence of two different universality 
classes. In order to test this result we study the force fluc- 
tuations of a disordered chain. The equilibrium length 
of the springs connecting the beads is chosen randomly 
(Poisson). The chain is then driven at constant velocity 
and the distribution of the slip sizes is calculated. The 
result is shown in Fig. ^ where we also report the distri- 
bution obtained with a periodic chain using similar pa- 
rameters. The two distributions are clearly indistinguish- 
able, casting some doubt on the relevance of disorder in 
the spring lengths. From this study, one would conclude 
that the two dynamics are in the same universality class. 



B. Constant force. 

For constant force driving, we employ system sizes 
varying from L = 20 to L — 200 and the density of pin- 
ning sites was chosen equal to unity [L — N) . In order to 
determine the exponent /3, we need an accurate estimate 
of the critical force, since an error in Fc can strongly bias 
the fit. 

Fig. H shows the value of the average velocity of the 
interface as a function of L. For F smaller than Fc in the 
limit of a large system v tends to zero, and for F greater 
than Fc and for the same limit {L ^ 1) v should tend a 
non vanishing value. In this way we can locate the crit- 
ical force, which results to be Fc = 2.195 ± 0.005. This 
result appears clearly from Fig. ||, from the log-log plot 
of u as a function of L. We see that F — 2.195 is compat- 
ible with a power-law behavior, whereas for F — 2.200 
in the large L limit the mean velocity tends to a nonzero 
constant. For F = 2.190 the velocity tends to zero faster 
than a power law in the limit of large system sizes. The 
numerical results are averaged over a number of disorder 
configurations which varies from 4000 for i = 20 to 500 
for L = 200. The curve plotted in Fig. p| allows also to 
estimate P/v = 0.16 ± 0.02 (see Eq. (ph)^ 



Next, we calculate the exponent f3 directly, plotting v 
vs (F - Fc). The fit in Fig. | yields (3 = 0.22 ± 0.02. 
We restrict the fit to the the six smaller values since a 
crossover to linear behavior is expected at high forces (i.e. 
V — F/rj) and this can bias the numerical estimate of the 
exponent. The simulations are made with a system of 
L = 540 particles, and the results are averaged over 100 
configurations of the disorder. In this way, we can obtain 
f3 and v as summarized in table I. 

To further test the consistency of our results, we cal- 
culate the exponent C measuring the scaling of W with L 
a.t F — Fc (see Eq. |lj) The results shown in Fig. |l^ give 
C = 1.28 ± 0.03 in agreement with the result obtained in 



Sec. IV A 



V. CONCLUSION 

In conclusion, we have investigated the dynamics of 
an elastic chain sliding on a disordered substrate. We 
have analyzed numerically the scaling close to the depin- 
ning transition focusing on the effect of different driving 
modes. Usually the problem is analyzed under constant 
force driving, while friction experiments are usually per- 
formed controlling the velocity. The two problems are 
closely related as discussed in Refs |,|,0|l9lj2|. We 
compute the critical exponents characterizing the tran- 
sition and analyze the effect of the driving velocity and 
the loading spring stiffness. Our results are in qualita- 
tive agreement with friction experiments performed with 
macroscopic asperities coupled by an elastic matrix. In 
general friction experiments our model will not apply 
since in many instances the Larkin length is extremely 
large . In addition, inertial effect are present in most 
cases and could lead to different force fluctuations. 

The exponents we measure agree well with the val- 
ues expected for the depinning of elastic interfaces in 
quenched disordered media [ p6| and with the renormal- 
ization group calculation of Ref. This is not sur- 
prising, since it is possible to show using the method dis- 
cussed in Ref. |lj] that the continuum limit of the model 
we study is described by 



dh{x, t) 
dt 



= DVh + F + fp{x,h), 



(14) 



where h{x,t) is a coarse grained version of Ui{t) and fp 
is a coarse grained random pinning force. It is thus ex- 
pected that the simulations performed in Ref. Q for a 
periodic chain agree with this result. The only difference 
between those simulations and ours lies in the way disor- 
der is implemented: in Ref. the pinning point are ar- 
ranged in a periodic structure and have random strength, 
while we use constant strength and random positions. 
We also note that the roughness exponent measured in 
Ref. ||] is very close to ( — 3/2, which is expected below 
the Larkin length, although the parameter employed do 
not seem to be consistent with that regime. 
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FIG. 2. Log-log plot of the size distribution of slip events 

for different values of V. In the inset we report the expo- 
nent r, obtained fitting the distributions, as a function of the 
driving velocity. 
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FIG. 4. Distribution of the slips events as a function of 
K showing that the cutoff increases when K decreases. The 
velocity V is kept constant. 
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FIG. 6. Width of the beads displacements as a function of 
K. The power law fit of the curve and the scaling relations 
imply that C = 1-26 ± 0.03. 
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FIG. 8. The variation of the mean velocity v of the chain 
as a function of the number of beads L and for various forces 
F. The inset shows the data for F far from the critical force 
Fc, and the main plot shows the data for F near the critical 
force, Fc ~ 2.195. 
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FIG. 9. Average velocity v as a function of the force F. The 
inset shows the v-F diagram. The sohd curve represents the 
result without pinning centers and the dotted line is the re- 
sult for the disordered substrate. The main graph represents 
the log-log plot of u as a function of _F — _Fc The exponent 
(5 is calculated using the six first points of the curve and the 
exponent is /? = 0.22 ± 0.02. 




FIG. 10. Width of the beads displacements W asa, function 
of L. The scaling relation between the two quantities gives 
C. We recover a value in good agreement with our previous 
results. 
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